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Percolation within Percolation 


Kit Tyabandha, Ph.D. 


The following programme was written in a macro language of Mat- 
lab. It was used in works leading to a Ph.D. thesis (Tiyapan, 2004). 
Percolation within percolation is the theme of that thesis. 


Regular 3-dimension lattice 


4 trr.m, regular 3-d tesselation, (c) Kit Tiyapan, 16 Dec 2002 
qn=size(q,1); 

in2n=size(in2,1); 

in3n=size(in3,1); 

ind5n=size(in5,1); 


p=; 

p{i,1}=[o0;0]’; 

p{1i,2}=sz; 

r=; 

s=(; 

t=O; 

for i=1:sz, 
r(1,i)=dx(m(i)) ; 
s(1,i)=dy(n(i)); 
t(1,i)=dz(z(i)); 

end 

p{i,3}=r’; 

p{i,4}=s’; 

9 p{1,5}=t’; 

20 iin=size(ii,1); 

21 tmp=ones(iin,1); 

22 tma=sparse(tmp,ii(:,1),tmp,1,sz); 

23 tape 

24 tmb=[] ; 

25 tmc=[]; 

26 tmd=[]; 

27 cnt=0; 

28 for i=1:sz, 

29 «if (~tma(i)) 


COCNDUFWNFOTOOANADAUTFLEWNEH 


30 cnt=cnt+1; 

31 tmp=([tmp;i,cnt] ; 
32 tmb=[tmb;r(i)]; 
33 tmc=[tmc;s(i)]; 
34 tmd=[tmd;t(i)]; 
38 end 

36 end 


37 p{2,1}=tmp; 

38 p{2,2}=cnt; 

39 p{2,3}=tmb; 

40 p{2,4}=tmc; 

41 p{2,5}=tmd; 

42 iiin=size(iii,1); 

43 tmp=ones(iiin,1); 

44 tma=sparse(tmp,iii(:,1),tmp,1,sz); 
45 tmp=[]; 
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tmb=([] ; 
tmc=[]; 
tmd=[]; 
cnt=0; 
for i=1:sz, 
if (~tma(i)) 
cnt=cntt1; 
tmp=([tmp;i,cnt] ; 
tmb=([tmb;r(i)]; 
tmc=[tmc;s(i)]; 
tmd=[tmd;t(i)]; 
end 
end 
p{3,1}=tmp; 
p{3,2}=cnt ; 
p{3,3}=tmb; 
p{3,4}=tmc; 
p{3,5}=tmd; 
iv=[ii,2*ones(iin,1)]; 
for i=1:iiin, 
tmp=0; 
for j=1:iin, 
if (“(iiiGi,1)-ii(j,1))) 
tmp=1; 
end 
end 
if (“tmp) 
iv=[iv;iii(i,:),3]; 
end 
end 
ivn=size(iv,1); 
tmp=ones(ivn,1); 
tma=sparse(tmp,iv(:,1),tmp,1,sz); 
tmp= al ; 
tmb=[] ; 
tmc=[] ; 
tmd=(] ; 
cnt=0; 
for i=1i:sz, 
if (~tma(i)) 
cent=cnt+1; 
tmp=[tmp;i,cnt] ; 
tmb=[tmb;r(i)]; 
tmc=[tmc;s(i)]; 
tmd=[tmd;t(i)]; 
end 
end 
p{4,1i}=tmp; 
p{4,2}=cnt; 
p{4,3}=tmb; 
p{4,4}=tmc; 
p{4,5}=tmd; 
vn=size(v,1); 
tmp=ones(vn, 1); 
tma=sparse(tmp,v(:,1),tmp,1,sz) ; 
tmp= al ; 
tmb=[] ; 
tmc=[]; 
tmd=[]; 
cnt=0; 
for i=1:sz, 
if (~tma(i)) 
cnt=cntt1; 
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tmp=([tmp;i,cnt] ; 
tmb=[tmb;r(i)]; 
tmc=[tmc;s(i)]; 
tmd=([tmd;t(i)]; 
end 
end 
{5,1}=tmp; 
Eis Obecne! 
p{5,3}=tmb; 
p{5,4}=tmc; 
p{5,5}=tmd; 
vi=[ii,2*ones(iin,1)]; 
for i=i:vn, 
tmp=0; 
for j=1l:iin, 
if (“(v(i,1)-ii(j,1))) 
tmp=1; 
end 
end 
if (~tmp) 
vi=[vi;v(i,:),5]; 
end 
end 
vin=size(vi,1); 
tmp=ones(vin,1); 
tma=sparse(tmp,vi(:,1),tmp,1,sz); 
tmp= ale 
tmb=[] ; 
tmc=(]; 
tmd=[] ; 
cnt=0; 
for i=1:sz, 
if (~tma(i)) 
cnt=cnt+1; 
tmp=([tmp;i,cnt] ; 
tmb=[tmb;r(i)]; 
tmc=[tme;s(i)]; 
tmd=[tmd;t(i)]; 
end 
end 
p{6,1}=tmp; 
p{6,2}=cnt; 
p{6,3}=tmb; 
p{6,4}=tmc; 
p{6,5}=tmd; 
vii=[iii,3*ones(iiin,1)]; 
for i=i:vn, 
tmp=0; 
for j=1:iiin, 
if (~(vGi,1)-iii(j,1))) 
tmp=1; 
end 
end 
if (~tmp) 
vii=[vii;v(i,:),5]; 
end 
end 
viin=size(vii,1); 
tmp=ones(viin, 1); 
tma=sparse(tmp,vii(:,1),tmp,1,sz); 
tmp= ae 
tmb=[] ; 
tmc=[]; 
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tmd=(] ; 
cnt=0; 
for i=1:sz, 
if (~“tma(i)) 
cnt=cntt1; 
tmp=(tmp;i,cnt] ; 
tmb=([tmb;r(i)]; 
tmc=[tmc;s(i)]; 
tmd=[tmd;t(i)]; 
end 
end 
p{7,1}=tmp; 
p{7,2}=cnt; 
p{7,3}=tmb; 
p{7,4}=tmc; 
p{7,5}=tmd; 
viii=iv; 
for i=1i:vn, 
tmp=0; 
for j=1:ivn, 
if (“(v(i,1)-iv(j,1))) 
tmp=1; 
end 
end 
if (“tmp) 
viii=([viii;v(i,:),5]; 
end 
end 
viiin=size(viii,1); 
tmp=ones(viiin, 1); 
tma=sparse(tmp,viii(:,1),tmp,i1,sz); 
tmp=[] ; 
tmb=[] ; 
tmc=(]; 
tmd=[] ; 
cnt=0; 
for i=1:sz, 
if (“tma(i)) 
cent=cntt1; 
tmp=(tmp;i,cnt] ; 
tmb=[tmb;r(i)]; 
tmc=[tmce;s(i)]; 
tmd=[tmd;t(i)]; 
end 
end 
p{8,1}=tmp; 
p{8,2}=cnt; 
p{8,3}=tmb; 
p{8,4}=tmc; 
p{8,5}=tmd; 
map=[]; 
tmp=ones(sz,1); 
map{1,1,1}=sparse(p{1,1}(:,1) ,tmp,p{1,1}(:,2),sz,1); 
cnt=p{1,2}; 


e=q; 
Ve [p{1,3} p14}, pC, 93] ; 
tmp=ones (p{2,2},1); 

for i=2:nx, 


map{i,1,1}=sparse(p{2,1}(:,1) ,tmp,cnt*#tmptp{2,1}(:,2),sz,1); 


for k=1:iin, 

map{i,1,1}(ii(k,1) ,1)=map{(i-1) ,1,1}(ii(k,2),1); 
end 
cnt=cnt+p{2, 2}; 
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235 % +iin; 
236 «©6©for k=1:qn, 
237 e=[e;map{i,1,1}(q(k,1)) ,map{i,1,1}(q(k,2))]; 


238 end 

239 for k=1:in2n, 

240 e=[e;map{i,1,1}(in2(k,1)),map{(i-1) ,1,1}(in2(k,2))]; 
24 end 

24 V=[V; (1-1) *dim1+tmp+p{2,3},p{2,4},p{2,5}]; 

243 end 


1 
2 
3 
4 tmp=ones(p{3,2},1); 
245 for j=2:ny, 
6 map{1,j,1}=sparse(p{3,1}(:,1) ,tmp,cnt*tmptp{3,1}(:,2),sz,1); 
7 for k=1:iiin, 
8 map{1,j,1}(iii(k,1),1)=map{1, (j-1) ,1} (iii(k,2) ,1); 
249 end 
250 cnt=cnt+p{3, 2}; 
251 % +iiin; 
252 for k=1:qn, 


253 e=[e;map{1,j,1}(q(k,1)) ,map{1,j,1}(q(k,2))]; 

254 end 

255 for k=1:in3n, 

256 e=[e;map{1,j,1}(in3(k,1)),map{1, (j-1) ,1} (in3(k,2))]; 
257 end 

258 «V=(V;p{3,3}, (j-1)*dim2*tmptp{3, 4} ,p{3,5}]; 

259 end 


260 tmp=ones(p{4,2},1); 

261 for i=2:nx, 

262 for j=2:ny, 

263 map{i,j,1}=sparse(p{4,1}(:,1),tmp,cnt*tmptp{4,1}(:,2),sz,1); 
264 for k=1:ivn, 


265 if (iv(k, 3) ==2) 

266 map{i,j,1}(iv(k,1) ,1)=map{(i-1) ,j,1} (ivék,2) ,1); 
267 else 

268 map{i,j,1}(iv(k,1) ,1)=map{i, (j-1) ,1} (iv(k,2) ,1); 
269 end 

270 end 

271 cnt=cnt+p{4,2}; 

272 % +ivn; 

273 for k=1:qn, 

274 e=[e;map{i,j,1}(q(k,1)) ,map{i,j,1}(q(k,2))]; 

275 end 

276 for k=1:in2n, 

277 e=[e;map{i,j,1}(in2(k,1)) ,map{(i-1),j,1} (in2(k,2))]; 
278 end 

279 for k=1:in3n, 

280 e=[e;map{i,j,1}(in3(k,1)) ,map{i, (j-1) ,1} (in3(k,2))]; 
281 end 


282 v=[V; (i-1) *dim1*tmp+p{4 ,3}, (j-1) *dim2*tmpt+p{4, 4}, p{4,5}]; 
283 end 

284 end 

285 tmp=ones(p{5,2},1); 

286 for i=2:nz, 

287 map{1,1,i}=sparse(p{5,1}(:,1) ,tmp,cnt*tmptp{5,1}(:,2),sz,1); 
288 for k=1:vn, 

289 map{i1,1,i}(v(k,1),1)=map{1,1, (i-1)} (v(k,2) ,1); 

290 end 

291 ent=cnttp{5,2}; 

292 for k=1:qn, 


293 e=[e;map{1,1,i}(q(k,1)) ,map{1,1,i} (q(k,2))]; 

294 end 

295 for k=1:in5n, 

296 e=[e;map{1,1,i}(in5(k,1)),map{1,1, (i-1)}(in5(k,2))]; 
297 end 
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298 =: V=([V;p{5,3},p{5,4}, (i-1) *dim3*tmptp{5, 53] ; 

299 end 

300 tmp=ones(p{6,2},1); 

301 for i=2:nx, 

302. for j=2:nz, 

303 map{i,1,j}=sparse(p{6,1}(:,1) ,tmp, cnt*tmpt+p{6,1}(:,2),sz,1); 
304 for k=1:vin, 


305 if (~ (wi (k,3)-2)) 

306 map{i,1,j}(vitk,1) ,1)=map{(i-1),1,j}(vitk,2) ,1); 
307 else 

308 map{i,1,j}(vi(k,1) ,1)=map{i,1, (j-1)} (vitk,2) ,1); 
309 end 

310 end 

311 cnt=cnt+p{6,2}; 

312 for k=1:qn, 

313 e=[e;map{i,1,j}(q(k,1)) ,map{i,1,j}(q(k,2))]; 

314 end 

315 for k=1:in2n, 

316 e=[e;map{i,1,j}(in2(k,1)) ,map{(i-1),1,j} (in2(k,2))]; 
317 end 

318 for k=1:in5n, 

319 e=[e;map{i,1,j}(ind(k,1)) ,map{i,1, (j-1)} (ind (k,2))]; 
320 end 

321 V=[V; (i-1) *dim1*tmp+p{6,3},p{6,4}, (j-1) *dim3*tmp+p{6,5}] ; 
322 end 

323 end 


324 tmp=ones(p{7 ,2},1); 

325 for i=2:ny, 

326 for j=2:nz, 

327 map{1,i,j}=sparse(p{7,1}(:,1),tmp,cnt*tmptp{7,1}(:,2),sz,1); 
328 for k=1:viin, 


329 if (~ (vii(k,3)-3)) 

330 map{1,i,j}(vii(k,1),1)=map{1, (i-1) ,j} (viitk,2),1); 
331 else 

332 map{1,i,j}(vii(k,1),1)=map{1,i, (j-1)} (viik,2),1); 
333 end 

334 end 

335 cnt=cnt+p{7 ,2}; 

336 for k=1:qn, 

337 e=[e;map{1,i,j}(q(k,1)) ,map{1,i,j}(q(k,2))]; 

338 end 

339 for k=1:in3n, 

340 e=[e;map{1,i,j}(in3(k,1)) ,map{1, (i-1) , j} (in3(k,2))]; 
341 end 

342 for k=1:in5n, 

343 e=[e;map{1,i,j}(ind5(k,1)) ,map{1,i, (j-1)} (ind (k,2))]; 
344 end 

345 v=[V;p{7,3}, (i-1) #dim2*tmp+p{7 ,4}, (j-1) *dim3*tmp+p{7 ,5}1; 
346 end 

347 end 

348 tmp=ones(p{8,2},1); 

349 for i=2:nx, 


350 for j=2:ny, 


351 for k=2:nz, 

352 map{i,j,k}=sparse(p{8,1}(:,1) ,tmp, cnt*tmptp{8,1}(:,2),sz,1); 
353 for m=1:viiin, 

354 if (~ (viii(m,3)-2)) 

355 map{i,j,k}(viii(m,1) ,1)=map{(i-1) ,j,k}(viii(m,2),1); 

356 elseif (~ (viii (m, 3) -3)) 

357 map{i,j.k}(viii(m,1) ,1)=mapf{i, (j-1) ,k} (viii(m,2),1); 

358 else 

359 map{i,j,k}(viii(m,1) ,1)=map{i,j, (k-1)}(viii(m,2) ,1); 

360 end 
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end 
cnt=cnt+p{8, 2}; 
for m=1:qn, 
entopmepttn) SiC ay »map{i,j,k}(q(m,2))]; 
en 
for m=1:in2n, 
e=[e;map{i,j,k}(in2(m,1)) ,map{(i-1) ,j,k}(in2(m,2))]; 
end 
for m=1:in3n, 
e=[e;map{i,j,k}(in3(m,1)) ,map{i, (j-1) ,k} (in3 (m, 2))]; 
end 
for m=1:in5n, 
e=[e;map{i,j,k}(in5(m,1)) ,map{i,j, (k-1)} (ind (m, 2))]; 
end 
V=[V; (i-1) *dim1*tmp+p{8, 3}, (j-1) *dim2*tmptp{8,4},... 
(k-1) *dim3*tmptp{8,5}] ; 
end 
end 
end 
en=size(e,1); 
Vn=size(V,1); 
figure(1); 
elt; 
hold on; 
for i=1:en, 


plot3(([V(e(i,1),1) ,V(e(i, 2) ,1)], [V(e(i, 1) ,2) ,V(e(i,2) ,2)],... 


[V(e(i,1),3) ,V(e(i,2) ,3)]); 
end 
axis off; 
axis equal; 
clf; 
hold on; 
tms=sum(nemat, 2); 
tmp=[] ; 
tma=17.2; 
% 5x5x5 units 
for i=1:en, 
if (tms (i) <tma) 
tmp=[tmp;e(i,:)]; 
end 
end 
tmn=size(tmp,1); 
for i=1:tmn, 
plot3(([V(tmp(i,1),1),V(tmp(i,2) ,1)] 


[v(tmp(i,1) ,2) ,V(tmp(i,2) ,2)], [v(tmp(i,1) ,3) ,v(tmp(i,2) ,3)]); 


end 
axis off; 
axis equal; 
4 for vertices 
nvmat=sparse (Vn, Vn) ; 
for i=1:en, 

nvmat (e(i,1),e(i,2))=1; 
nvmat (e(i,2) ,e(i,1))=1; 
end 
A=V; 
N=Vn ; 
lmat=sparse(1,N); 
umat=sparse(1,N) ; 
LB=min(V(:,1)); 
UB=max(V(:,1)); 
rng=UB-LB; 
LBv=.05*rng+LB ; 
UBv=UB-LBv ; 
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for i=1:Vn, 
if (V(i, 1)<LBv) 
lmat (1,i)=1; 
end 
if (V(i,1)>UBv) 
umat (1,i)=1; 
end 
end 
nmat=nvmat ; 
Blocked=randperm(Vn) ; 
[pc, cord, tsries]=perc(N, 1mat,umat ,nmat) ; 
4 for edges 
evm=sparse (en, Vn) ; 
for i=1:en, 
evm(i,e(i,1))=1; 
evm(i,e(i,2))=1; 
end 
nemat=sparse (en, en) ; 
for i=1:Vn, 
tmp=find(evm(: ,i)); 
tmn=size(tmp,1); 
for j=1:(tmn-1), 
for k=(j+1):tm, 


nemat (tmp(j) ,tmp(k))=1; 
nemat (tmp(k) ,tmp(j))=1; 
end 
end 
end 
A=e; 
N=en; 


lmat=sparse(1,N); 
umat=sparse(1,N) ; 
for i=1:N, 
if ((V(ACi, 1) ,1)<=LBv) | (V(A(i,2) ,1)<=LBv)) 
Imat (1,i)=1; 
elseif ((V(AC(i,1),1)>=UBv) | (V(A(i,2) ,1)>=UBv)) 
umat (1,i)=1; 
end 
end 
nmat=nemat ; 
Blocked=randperm(N) ; 
[pc, cord, tsries]=perc(N,1mat, umat ,nmat) ; 
4 cep 
clear all; 
sz=9; 
nx=5; 
ny=5; 
nz=5; 
r=1; 
tmp=sqrt (2) +r; 
dx=[0,2*r,4*r]; 
dy=[0,2*r,4+#r]; 
dz=(0,tmp,2*tmp] ; 
dimi=max (dx) ; 
dim2=max (dy) ; 
dim3=max (dz) ; 


31,531,6;2,533,534,5;5,63;5,7;5,8;5,9]; 
m=[1,3,1,3,2,1,3,1,3]; 
n=[1,1,3,3,2,1,1,3,3]; 
z=[1,1,1,1,2,3,3,3,3];0=[1,2,3,4,5,6,7,8,9]; 
ii=[1,2;3,4;6,7;8,9]; 
iii=[1,3;2,4;6,8;7,9]; 
v=[1,6;2,7;3,8;4,9]; 
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487 in2=[]; 
488 in3=[]; 
489 ind=(]; 


Bibliography 


K N Tiyapan. Division of space by Voronoi graphs, percolation within 
percolation and application to the models of porous membranes. PhD 
Thesis. University of Manchester, 2004 


Kit Tyabandha. Percolation within Percolation, Theory. Tyabandha Journal 
of Arts and Science. 3, 2, 2006 


270 July 2006 Vaen Sryayudhya, Editor 


